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We study front propagation and diffusion in the reaction-diffusion system A ^ A + A on a 
lattice. On each lattice site at most one A particle is allowed at any time. In this paper, we analyze 
the problem in the full range of parameter space, keeping the discrete nature of the lattice and 
the particles intact. Our analysis of the stochastic dynamics of the foremost occupied lattice site 
yields simple expressions for the front speed and the front diffusion coefficient which are in excellent 
agreement with simulation results. 
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INTRODUCTION 

In this paper, we study the propagation and diffusion 
of a front in the A ^ A + A reaction on a chain, in the 
case that there cannot be more than one A-particle on 
each lattice site ( "hard-core exclusion" ) . The front prop- 
agation problem we consider is the following. We start 
from a situation illustrated in Fig. ^a) in which there 
are no A-particlcs at all on the right half of the system, 
while there is a nonzero density of particles on the left. 
The object of study is then the asymptotic average speed 
V with which the region with a nonzero density of parti- 
cles expands to the right, as well as the effective diffusion 
coefficient -D/ of this "front" . For the hard-core exclusion 
problem, the front position is most conveniently defined 
as the position of the foremost (rightmost) particle, see 
Fig. ^a-6j. The average front speed and front diffusion 
coefficient are then the average drift speed v and the dif- 
fusive spreading \/Dft of the width of the probability 
distribution Pkj- (t) for the location % of the foremost oc- 
cupied lattice site, as illustrated in Fig. Hfcj. One of the 
main results of the paper is a simple expression for v and 
Df, which is accurate in the range where the deviations 
from the mean field theory are large. Our results reduce 
to an exact expressions derived before for the particular 
case in which the particle diffusion coefficient D and an- 
nihilation rate W are equal 1\ and our expression for the 
front speed v reduces to the approximate expression ob- 
tained for the special case = in j3, 0; 01 ■ addition, 
we study the average particle profile behind the foremost 
occupied lattice site and analyze how its behaviour af- 
fects the average front speed and diffusion. 

The perspective of this work lies in the issues that have 
emerged from the surprising findings for fronts in this 
reaction-diffusion system in the limit in which N , the av- 
erage number of particles per lattice site in equilibrium, 
is large. In a lattice model, one can tune N by allowing 
more than one particle per site (no hard-core exclusion) 
and changing the ratio fcb/fcd, where fcb is the reaction 
rate for birth processes A 2A and fed the reaction rate 
for death processes, 2A A, as the average equilibrium 
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FIG. 1: (a) The type of initial condition we consider for our 
stochastic model, (b) Illustration of a typical snapshot of the 
state of the system at finite time. The foremost particle has 
advanced to the right relative to the one where it started at 
t = 0. (c) Qualitative sketch of the probability distribution 
function for the foremost particle at t = (dashed line) at 
large times t; the center of the peak drifts with speed v, while 
the peak widens proportional to \/Df t. 



number of particles N = fcb/fcd. In the limit iV— >oo, 
the normalized particle density pi = Ni/N then obeys a 
mean- field equation which is a lattice analogue of the con- 
tinuum reaction diffusion equation dtp — Dd^p + p — p^ , 
where D is the diffusion rate of individual particles on 
the chain. The front problem mentioned above, i.e., the 
propagation of a front into an empty region, then corre- 
sponds in the mean-field limit N—^oo to a front propa- 
gating into the linearly unstable state p = (the mean 
field behavior is also obtained in the limit in which the 
particle diffusion coefficient D oo |^ H, 0| , but we will 
focus on the case in which the diffusion coefficient is finite 
and comparable to the growth and annihilation rates). 
The behavior of such fronts in deterministic continuum 
equations has been studied since long and is very well 
understood (see, e.g., 0,0)- Since the nonlinear front 
solutions are essentially "pulled along" by the growth of 
the leading edge where p <C 1, such fronts are often re- 
ferred to as pulled fronts The remarkable discovery 
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of the last few years has been that since the propagation 
is driven by the region where p is small, they are par- 
ticularly sensitive to the discrete nature of the particles 
which manifests itself in changes in the dynamics when 
p becomes of order Indeed, Brunet and Derrida 

discovered that the convergence to the mean-field limit 
is extremely slow with N: the average front speed v con- 
verges as 1/ In^ N to the mean- field value (21] • This is in 
contrast to the fact that for pushed fronts, the conver- 
gence to asymptotic speed behaves as a power of \/N. 
This slow convergence has been confirmed for a variety 
of models In addition^ in a 

model that Brunet and Derrida studied in Ref. the 
front diffusion coefficient Df was numerically shown to 
vanish only as 1/ In^ N . 

The dominant asymptotic correction to the mean-field 
result for the front speed in the limit N oo traces 
simply to the change in the dynamics at p = 0{1/N) 
0, and as a result appear to be universal. However, 
all corrections beyond the asymptotic one appear to de- 
pend non-universally on the detailed stochastic dynam- 
ics at the foremost occupied site and those closely behind 
it, where asymptotic techniques are of no use since the 
number of particles involved in the dynamics is small ■ 
Moreover, the stochastic dynamics in the tip region even 
seems to be strongly nonlinearly coupled to the uniformly 
translating average front profile behind the tip. 

For analyzing these effects for finite values of the par- 
ticle diffusion coefficient D and particle number N, it is 
found to be expedient to develop a stochastic front de- 
scription by focussing on the behavior of the foremost 
particle or the foremost occupied bin 0|. As it turns 
out, this idea traces back to the earlier work by Kerstein 
0, 3 and Bramson and coworkers Q . These authors an- 
alyzed the average front speed v for a special case of the 
model we investigate here, namely the case in which the 
particle annihilation rate W = 0. In this case, one can 
formulate a self-consistent dynamics for the two foremost 
particles |l5|. but this important simplification is lost 
when W ^ . Motivated by the desire to understand 
the ingredients necessary to analyze the stochastic front 
behavior for finite values oi D, W and N, we focus here 
on analyzing both v and Df in the case in which all the 
transition rates are comparable; our analysis includes the 
special point D = W where an exact result was obtained 
by ben-Avraham 0. 



THE MODEL, FRONT SPEED AND FRONT 
DIFFUSION 



We now turn to the details of our model and our re- 
sults for the stochastic fronts. We consider a chain on 
which A particles can undergo the following three basic 
moves, shown in Fig. [3 ('j^ A particle can diffuse to any 
one of its neighbour lattice sites with a diffusion rate D, 
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FIG. 2: The microscopic processes that take place inside the 
system: (i) A diffusive hop with rate _D to a neighboring 
empty site; (ii) Creation of a new particle on a site neigh- 
boring an occupied site with rate e; (iii) Annihilation of a 
particle on a site adjacent to an occupied site at a rate W. 



provided this neighboring site is empty, (ii) Any particle 
can give birth to another one on any one of its empty 
neighbour lattice site with a birth rate e. (iii) Any one 
of two A particles belonging to two neighbouring filled 
lattice sites can get annihilated with a death rate W. 

Note that in the above formulation, diffusive hops to 
neighboring sites which are occupied are not allowed. We 
can also think about these stochastic moves differently: 
for example, we can allow nearest neighbour diffusive 
hops to a site which is already occupied be followed by 
an instantaneous annihilation of one of the two particles. 
If we do so, then the diffusive process contributes to the 
annihilation of particles. However, in this paper we shall 
stick to the convention that diffusive hops are allowed 
only to empty sites. 

As noted before, earlier work on models of this type 
includes that of Kerstein 0, 0| and Bramson et al. Q 
on the case W = and that of ben-Avraham on the 
case Z3 = 0] (also, variants of this model have been 
analyzed in Refs. [13; llM llSl ) ■ Notice that in the general 
case there are essentially only two nontrivial parameters 
in the model, e.g. the ratios D/e and D/W, since an 
overall multiplicative factor just sets the time scale. Our 
interest is in the parameter range where both of these 
ratios are 0(1); when these ratios tend to infinity, the 
front speed approaches the mean field value 0, 0j ■ 

For an ensemble of front realizations, let us denote the 
probability distribution for the foremost occupied lattice 
site to be at lattice site kf hy Pkj{t). The evolution of 
Pkf (t) is then described by 



dt 



^ = [D 
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Here P^^'^it) and Pfc"P*^(t) denote the joint probabili- 
ties that the foremost particle is at site kf and that the 
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site fc/ — 1 is occupied or empty, respectively. Clearly, 
^fc, it) = Pk^it) + PkT'it)^ Ek, Pk, it) = 1. The 



first term on the r.h.s. of Eq. describes the increase 
in (t) due to the advancement of a foremost occupied 
lattice site from position A;/ — 1, while the second term 
describes the increase in Pkfit) due to the retreat of a 
foremost occupied lattice site from position kj+l. The 
third and the fourth terms, respectively, describe the de- 
crease in Pkj (t) due to the advancement and retreat of a 
foremost occupied lattice site from position kf. 

From the definition of P^f (t) , the mean position and 
the width of the distribution for the positions of the 
foremost occupied lattice sites are defined as x{t 



E,, kfP,^{t) and {AxHt))^^j:,^[kf-xit)]^PkAt) 
The mean speed and diffusion coefRcient of the front are 
thus given in terms of these quantities as the i— >oo limit 
ofv = dx{t)/dt and (Aa;2(t)) = 2D ft — see Fig. ^c). 
To obtain them, we need the expressions of P^^'^ (t) and 
P™''*-*'(i). To start with, we have 



P^-(i) = Pfe,_iP,,(i), 



(2) 



where Pkf-i is the conditional probability of having the 
(/c/ — l)th lattice site occupied (the foremost particle is 
at the kfth lattice site). The set of conditional occu- 
pation densities pkf~m for m > 1 can be thought of as 
determining the front profile in a frame moving with each 
front realization. For obtaining v and D f , we simply need 
to know the asymptotic long-time limit pkf-iit — > oo), 
which from here on we will denote simply as pkj - 1 . Given 
Pkj-i, it is then straightforward to obtain from Eq. 
and the conditions Pkj{t) = Pk'f'^it) + P^^^^^it) and 

dx 

V ^ — ^ e - pfe - D) and 

at ' 



d{Ax^) 
Jt 



2D + e + pkf-i{W - D) . (3) 



Of these, the second equation indicates that the front 
wandering is diffusive, and an expression of the front dif- 
fusion coefficient Df is therefore given by 

Df ^^[2D + e + pk,-iiW - D)] . (4) 

As noted already by ben-Avraham |jy|] in a continuum 
formulation of the present model, for the special case 
D = W the unknown quantity Pkf-i drops out of Eq. 
©; it thus leads to the exact results v = e and Df — 
D + e/2 a.s a. special cases of Eq. for = W. We also 
note that if we use Eq. (0) in Eq. the latter equation 
has the form of the master equation for a single random 
walker on a chain. Thus we can think of the foremost 
particle as executing a biased random walk, and Df as 
the effective diffusion coefficient of this walker. Moreover, 
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FIG. 3: Comparison of tlie expression of v in Eq. @ (solid 
line) with stochastic simulation data (filled circles), for D = 
e = 0.25. The error in the data is of the order of the size of 
the symbols. The corresponding data point for D — W , as 
analyzed in Ref. 0] is shown by the larger open circle. 



if we eliminate Pkf-i from Eqs. 
following exact relation 



© and Q), we get the 



v/2 + Df = D + e . 



(5) 



In order to obtain an explicit prediction for v and Df, 
we need an expression for pkj-i- Far behind the front 
the particle density will approach the homogeneous equi- 
librium density p: limm^oo Pkj-m = P- From the master 
equation it is easy to show that the homogeneous equi- 
librium solution for the total probability is of product 
form (so that the probability of having different sites is 
occupied is uncorrelated) , and that the equilibrium oc- 
cupation density p is simply given by p = e/{£ -\- W). 

The crudest approximation for the front profile pkj~m 
and in particular for Pkf-i is to just take Pkj-i ~ p. 
Substitution of this approximation into Eqs. and Q 
immediately yields our main result. 



e{e + D) 
e + W 



and Df = 



(e + 2W){D + e) 
2{s + W) 



(6) 



For W = 0, the expression for v reduces to the one ob- 
tained in [illH. 

In what follows, we will first compare these approxi- 
mate expressions for v and Df to the results of computer 
simulation for the case D / e = \ , and then investigate the 
appropriateness and shortcomings of the approximation 
Pki-i ~ p. 

The comparison of Eq. ((BJ with stochastic simulation 
data for D = e = 0.25 are presented in Figs. |21 and 
01 as a function of W ioY D = e = 0.25. The simula- 
tion algorithm has been adopted from and is es- 
sentially the same one as in 9]. The speed v has been 
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FIG. 4: Comparison of the front diffusion coefficient according 
to Eq. Q (solid line) with stochastic spreading data (filled 
circles) and with Eq. @ (open triangles), for D = e — 0.25. 
The large open circle once again corresponds to the direct 
measurement of the effective front diffusion coefficient for D = 
W, as analyzed in Ref. Q. 



obtained directly from the average position of the fore- 
most occupied lattice site in a single long run accord- 
ing to v{t) = [x{t) — x{to)]/{t — to) corresponding to 
x{t) — x{tQ) — 15000 consecutive forward jumps. The 
diffusion coefficient has been determined both from the 
speed measurements via ^ and from data for the aver- 
age diffusive spreading during 1000 time intervals At up 
to 500 taken from 5 long runs (of which the data from the 
first 5000 consecutive forward jumps of the foremost oc- 
cupied lattice site were ignored, so as to eliminate initial 
transient effects). For each of these runs, the mean square 
displacement (Ax^) was confirmed to grow linearly with 
time. Figures 13 and 01 show that our approximate ex- 
pressions (jnj for the speed and diffusion coefficient (solid 
line) are quite accurate for D/e = 1 over the whole range 
of values of W where we have performed simulations; 
interestingly, the values of Df obtained from the speed 
measurements via Eq. (j^l are more accurate than those 
obtained directly from the diffusive spreading. The error 
bars in Fig. ^correspond to the standard deviations of 
Df values obtained from 5 long runs. 

We now return to the issue of the appropriateness of 
the assumption pkf-i = P- While the agreement between 
the theoretical prediction for v and Df gives empirical 
evidence that this assumption is a reasonably good one, 
we see from Fig. |31 that although Eq. Q agrees well 
with the simulation data, there are small but systematic 
deviations on both sides of this region. These deviations 
can be explained as follows: As Vl^ ^ 0, p t 1- f^^i' 
behind the front all lattice sites are occupied. However, 
the density of particles just behind the foremost one is 



FIG. 5: Relative deviation d = {pk — p)/p of the average 
density from p ~ sjie 4- W) for the first six lattice sites to the 
left of the foremost occupied lattice site. A;/, for D = e = 0.25 
and three different values of W . 

smaller, since it takes a finite time for the density to relax 
to the asymptotic one. For large values of W, the effective 
diffusion rate is much larger than the drift rate, as Eq. 
© shows. As a result, once again the density of particles 
just behind the foremost one also has relatively small 
time to relax to the asymptotic value. This is reflected 
in the difference between Pkf-i and p in Fig. 

The above trends are borne out by the simulation re- 
sults of Fig. |S1 where we plot the relative deviation 
d = [pk — p)/p for k = kf ~ 1, . . . , kf — 6. First of all, 
the data confirm that unless the value W is too small, 
Pkj-i — J) is quite a good approximation, and that the 
density behind the foremost particle is enhanced for large 
W and reduced for small W. We also note that we have 
verified that if one substitutes the Pkf-i values for W = 
and W ~ 0.8 from Fig. (Sjinto Eq. one does recover 
the corresponding measured speeds, as one should. 



CONCLUSION 

In conclusion, this work clearly illustrates that the con- 
cept of the dynamics of the foremost occupied lattice site, 
in Refs. y, ^ Q| and here, can be a viable route to- 
wards analyzing the front propagation and diffusion in 
stochastic lattice models. In the present < 1 model a 
simple approximation for the interaction of the foremost 
particle with the front region behind it already yields 
quite accurate results for v and Df. We hope that this 
success provides new motivation and inspiration to tackle 
the complicated case in which N is large but finite. 

In principle, it should be possible to extend the anal- 
ysis in the spirit of the one developed by Kerstein 0, 
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to get successively more accurate expressions for Pkf-i, 
and correspondingly for the front speed and diffusion co- 
efficient. In particular, such extensions might allow one 
to use the results in a wider parameter range, such as 
D/W oo while D/e - or D/e oo while 

W/e ^ 0{1). However, inspection of the earlier analysis 
suggests that such higher order analytical expressions of 
Pkj-i are less trivial to obtain than one mi ght expect at 
first sight. More precisely, in the light of [l^, lla|> it is 
clear that for W ^ 0, the master equation for the proba- 
bility that the two foremost particles are separated by k 
lattice sites couples to probability distributions involving 
particles that are further back. While it is certainly pos- 
sible to solve the master equation numerically, it does not 
appear to lead one to an analytical expression of Pkf~i 
that provides a better approximation than what we have 
used in this paper. 

The work by D. P. and that of G.T. during an earlier 
stay at Universiteit Leiden was supported by the Foun- 
dation FOM ( "Fundamenteel Onderzoek der Materie"). 
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